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Summary. Infiation offers a simple model for very early evolution of our Universe 
and the origin of primordial perturbations on large scales. Over the last 25 years 
we have become familiar with the predictions of single-field models, but inflation 
with more than one light scalar field can alter preconceptions about the inflationary 
dynamics and our predictions for the primordial perturbations. I will discuss how 
future observational data could distinguish between inflation driven by one field, or 
many fields. As an example, I briefly review the curvaton as an alternative to the 
inflaton scenario for the origin of structure^ 



1 Introduction 

Inflation provides an attractively simple model for the early evolution of our Uni- 
verse, which can produce a large, spatially flat and largely homogeneous observable 
universe. It also provides a source for small primordial perturbations which are the 
origin of the large-scale structure in our Universe today. The vacuum fluctuations 
of any light scalar field present during inflation can be swept up by the inflationary 
expansion to scales much larger than the Hubble scale. 

Inflation is most commonly discussed in terms of a potential energy which is a 
function of a single, slowly rolling, scalar field. Single field, slow-roll inflation pro- 
duces an almost Gaussian distribution of adiabatic density perturbations on super- 
Hubble scales with an almost scale-invariant spectrum. But supersymmetric field 
theories can contain many scalar fields that could play a role during inflation and 
string theory, and other higher-dimensional theories, yield four- dimensional effec- 
tive actions with many moduli fields describing the higher dimensional degrees of 
freedom. One should be aware of the different possibilities that open up in particle 
physics models containing more than one light scalar field during inflation. 

The presence of multiple fields during inflation can lead to quite different infla- 
tionary dynamics, that might appear unnatural in a single field model, and to spectra 
of primordial perturbations that would actually be impossible in single field models. 



^ Based on a talk given at the 22nd lAP Colloquium, "Inflation -1-25", Paris, June 
2006 
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The presence of multiple light fields during infiation leads to the generation of non- 
adiabatic field perturbations during infiation. This can alter the evolution of the 
overall curvature perturbation, for instance leading to detectable non-Gaussianity, 
and may leave residual isocurvature fluctuations in the primordial density pertur- 
bation on large scales after infiation, which can be correlated with the curvature 
perturbation. Such alternative models are interesting not only as theoretical possi- 
bilities, but because they could be distinguished by increasingly precise observations 
in the near future. 

In this short review I will discuss some of the distinctive observational predictions 
of inflation in presence of more than one scalar field. For a more comprehensive 
review of infiationary dynamics and reheating with multiple fields see Ref . [1] . 



2 Homogeneous scalar field dynamics 

The time-evolution of a single, spatially homogeneous scalar field is governed by the 
Klein-Gordon equation 

4> + 3H(j} = -^, (1) 



where the Hubble expansion rate is given by the Friedmann constraint 
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(2) 



Multiple scalar fields obey the Klein-Gordon equation 

^j + 3Hipi = ^^iy"Uj] , (3) 




where one must allow for the possibility that the potential energy is given by a sum 
over many terms 

V = Y,Uj. (4) 

J 

The wider range of interaction potentials possible in multiple field models leads to 
possibilities such as hybrid inflation. 

In hybrid infiation models [2l|3] the inflaton field, ipi, can roll towards the non- 
zero minimum of its potential, Ui = Vo + m\ifi\/2, which would lead to eternal 
inflation into the future in a single field model. But in a hybrid model there is a 
second waterfall scalar field trapped during inflation in a local minimum, ip2 = 0, 
with a potential, e.g., U2 = {g'^ip\—m%)tp\/2 which becomes unstable below a critical 
value of the (p\ field, triggering an instability of the vacuum energy driving infiation, 
rapidly bringing inflation to an end. 

Another more subtle change enters through the Friedmann constraint: 



Even in the absence of explicit interactions in the scalar fleld Lagrangian, the flelds 
will still be coupled gravitationally. In particular the Hubble expansion rate that 
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enters the Klein-Gordon equation Q is due to the sum over all fields in Eq. (O and 
this can also alter the field dynamics even if the potential for each individual field is 
left unchanged. The additional Hubble damping present due to multiple fields can be 
used to drive slow-roll infiation in assisted infiation models [1] where the individual 
potentials would be too steep to drive infiation on their own. 

The original assisted infiation model [4] considered n scalar fields with steep 
exponential potentials 

V = Uio exp {~Xi(pi/mp) (6) 

where I have used the reduced Planck mass nip — (StvG)^^ . Each scalar field poten- 
tial is too steep to drive inflation on its own if Xj > 2, but the additional damping 
effect due to the presence of the other scalar fields leads to a particular power-law 
infiation solution, a oc f where p = where the combined fields have an effective 
potential V oc exp{—Xa/mp) with 

I ' 

Thus A for many fields as n —* oo and we can have slow-roll inflation even when 
each Af > 2. 

Even though the background dynamics can be reduced to an equivalent sin- 
gle field with a specified potential 5 , there is an important qualitative difference 
between the inflationary dynamics in multiple field infiation with respect to the sin- 
gle field case. The Hubble damping during inflation drives a single scalar fleld to 
a unique attractor solution during slow-roll inflation where the Hubble rate, fleld 
time- derivative and all local variables are a function of the local field value: H{(f>), 
4'{4') 1 etc. This means that the evolution rapidly becomes independent of the initial 
conditions. 

In multiple field models we may have a family of trajectories in phase space 
where, for example, the Hubble rate at a particular value of ipi is also dependent 
upon the value of ip2- In this case the infiationary dynamics, and hence observational 
predictions, may be dependent upon the trajectory in phase space and thus the 
initial field values. It is this that allows non-adiabatic perturbations to survive on 
super-Hubble scales in multiple field inflation. 

It is important to distinguish here between models, such as most hybrid models, 
with multiple flelds but only one light direction in field space with small effective 
mass d^V/d(p^ ^ during infiation, and models with many light fields, such as 
assisted inflation models. Only models with multiple light flelds can have multiple 
slow-roll trajectories. 



2.1 Inflaton field direction during inflation 

It is convenient to identify the inflaton field direction as the direction in field space 
corresponding to the evolution of the background (spatially homogeneous) fields 
during infiation [6] (see also [3 IB])- Thus for n scalar fields ipi, where I runs from 1 
to n, we have 

a = / ^ anpidt , (8) 
•' I 



4 David Wands 



where the inflaton direction is defined by 

The n evolution equations for the homogeneous scalar fields Q can then be 
written as the evolution for a single infiaton field ([1} 

a + 3H& + V^^0, (10) 

where the potential gradient in the direction of the infiaton is 

The total energy density and pressure are then given by the usual single field results 
for the infiaton. 



2.2 An example: Nflation 

A topical example of multiple field infiation is Nfiation. Dimopoulos et al proposed 
this model based on the very large number of axion fields predicted in low energy 
effective theories derived from string theory. Near the minimum of the effective 
potential the fields have a potential energy 

F=i^m?^?. (12) 

This form of potential with a large number of massive fields was also previously 
studied by Kanti and Olive [10] and Kaloper and Liddle [llj . 

With a single scalar field the quadratic potential yields the familiar chaotic 
infiation model with a massive field V = m^0^/2. But to obtain infiation with a 
single massive field the initial value of the scalar field must be several times the 
Planck mass and there is a worry that we have no control over corrections to the 
potential at super-Planckian values in the effective field theory 12 . But with many 
scalar fields the collective dynamics can yield infiation even for sub-Planckian values 
if there are a sufficiently large number of fields. Kim and Liddle [13) have found that 
for random initial conditions, —nip < tpi{0) < rup, the total number of e-folds is 
given by n/12, where n is the total number of fields. Thus we require n > 600 for 
sufficient infiation if none of the fields is allowed to exceed the Planck scale. This 
may seem to be a large number, but Dimopoulos et al cite string theory models 
with of order 10^ axion fields. 

As remarked earlier, in the presence of more than one light field, the trajectory 
in field space at late times, and hence the observable predictions, may be dependent 
upon the initial conditions for the different fields. But Kim and Liddle [TS] found 
evidence for what they called a "thermodynamic" regime where the predicted spec- 
tral index, niz, for the primordial curvature perturbations that arise from quantum 
fiuctuations of the scalar fields, became independent of the precise initial conditions 
for a sufficiently large number of fields. In fact infiation with an arbitrary number 
of massive fields always yields a robust prediction for the tensor-to-scalar ratio r in 
terms of the number of e-foldings. A'', from the end of infiation [12] 
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r = 8/iV, (13) 

completely independently of the initial conditions. Thus Nflation seems to be an 
example of a multiple field model of infiation which makes observable predictions 
which need not depend upon the specific trajectory in field space. 

In the limit where the masses become degenerate, mf — > m, the Nflation dy- 
namics becomes particularly simple. The fields evolve radially towards the origin 
and the potential (|12p reduces to that for a single field 

V ^ im^a^ , (14) 

where a is the inflaton field ([8]). Thus in this limit Nflation reproduces the single 
field prediction for the tensor-scalar ratio r = 0.16 and the spectral index n-n = 0.96. 
However the presence of n light fields during Nfiation also leads to n — 1 isocurvature 
modes during inflation and these have an exactly scale-invariant spectrum (up to 
flrst order in the slow-roll parameters) in the limit of degenerate masses [TJ]. In the 
following sections I will describe some of the distinctive predictions that can arise 
due to the existence of such non-adiabatic perturbations during infiation. 



3 Primordial perturbations from inflation 

I have so far only presented equations for the dynamics of homogeneous scalar 
fields driving inflation. But to test theoretical predictions against cosmological ob- 
servations we need to consider inhomogeneous perturbations. It is the primordial 
perturbations produced during inflation that offer the possibility of determining the 
physical processes that drove the dynamical evolution of the very early universe. In 
the standard Hot Big Bang model there seems to be no way to explain the existence 
of primordial perturbations, during the radiation dominated era, on scales much 
larger than the causal horizon, or equivalently the Hubble scale. But infiation takes 
perturbations on small, sub-Hubble, scales and can stretch them up to arbitrarily 
large scales. 

3.1 Scalar field perturbations, without interactions 

Consider an inhomogeneous perturbation, i^i fiit) + Sipi{t,x.), of the Klein- 
Gordon equation (|3]) for a non-interacting scalar field in an unperturbed FRW uni- 
verse. (I will include perturbations of the metric and other fields later, but for 
simplicity I will neglect this complication for the moment.) 

Sifj + 3H5ifj + (m? - V^) Sifii ^0, (15) 

where the effective mass-squared of the field is mj — d'^V/dip^, and is the 
spatial Laplacian. Decomposing an arbitrary field perturbation into eigenmodes of 
the spatial Laplacian (Fourier modes in flat space) V'^Sipi — —{k'^ /a'^)Sipi, where 
k is the comoving wavenumber and a the FRW scale factor, we find that small- 
scale fiuctuations in scalar fields on sub-Hubble scales (with comoving wavenumber 
k > aH) undergo under-damped oscillations, and on sufficiently small scales are 
essentially freely oscillating. Normalising the initial amplitude of these small-scale 
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fluctuations to the zero-point fluctuations of a free fleld in flat spacetime we have 

m 

— ikt I a 

av2fc 



During an accelerated expansion a — aH increases and modes that start on sub- 
Hubble scales {k > aH) axe stretched up to super-Hubble scales (k < aH). Perturba- 
tions in light fields (with effective mass-squared < 9H^ /A) become over-damped 
(or "frozen-in" ) and Eq. (|16|) evaluated when k ~ aH gives the power spectrum for 
scalar field fluctuations at "Hubble-exit" 

_ 47rfc^ i„ ,1 /H 



(27r) 



Heavy fields with > 9H^ /A remain under-damped and have essentially no per- 
turbations on super-Hubble scales. But light fields become over-damped and can 
be treated as essentially classical perturbations with a Gaussian distribution on 
super-Hubble scales. 

Thus infiation generates approximately scale-invariant perturbation spectra on 
super-Hubble scales in any light field (for < H'^ and < H'^). 



3.2 Scalar fleld and metric perturbations, with interactions 

The simplified discussion in the preceding subsection gives a good approximation to 
the scalar field perturbations generated around the time of Hubble-exit during slow- 
roll infiation, where field interactions and metric backreaction are small. However 
to accurately track the evolution of perturbations through to the end of infiation 
and into the radiation dominated era we need to include interactions between fields 
and, even in the absence of explicit interactions, we need to include gravitational 
backreaction. 

For an inhomogeneous matter distribution the Einstein equations imply that 
we must also consider inhomogeneous metric perturbations about the spatially fiat 
FRW metric. The perturbed FRW spacetime is described by the line element [341 [1] 

ds^ = -(1 + 2A)dt^ + 2adiBdx'dt 

W [(1 - 2V')(5ij + 2dijE + h,j] dx'dx^ , (18) 

where di denotes the spatial partial derivative d/dx^. We will use lower case latin 
indices to run over the 3 spatial coordinates. 

The metric perturbations have been split into scalar and tensor parts according 
to their transformation properties on the spatial hypersurfaces. The field equations 
for the scalar and tensor parts then decouple to linear order. Vector metric pertur- 
bations are automatically zero at first order if the matter content during infiation 
is described solely by scalar fields. 

The tensor perturbations, htj, are transverse {d^hij = 0) and trace-free {5'^hij — 
0). They are automatically independent of coordinate gauge transformations. These 
describe gravitational waves as they are the free part of the gravitational field and 
evolve independently of linear matter perturbations. 

We can decompose arbitrary tensor perturbations into eigenmodes of the spatial 
Laplacian, V^e^j = ~{k^ /a^)eij, with comoving wavenumber fc, and scalar ampli- 
tude h{t): 
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h.,=h{t)e^'-\x), (19) 

with two possible polarisation states, + and x . The Einstein equations yield a wave 
equation for the amplitude of the tensor metric perturbations 

h + 3Hh+ —h^O, (20) 

This is the same as the wave equation (|15|l for a massless scalar field in an unper- 
turbed FRW metric. 

The four scalar metric perturbations A, diB, ipSij and dijE are constructed from 
3-scalars, their derivatives, and the background spatial metric. The intrinsic Ricci 
scalar curvature of constant time hypersurfaces is given by 

(3)ii=±vV- (21) 
a-' 

Hence we refer to i/j as the curvature perturbation. 

First-order scalar field perturbations in a first-order perturbed FRW universe 
obey the wave equation [l] 

Sifj + 3H5(pj H — TrSipi + VijSifj 
a-' 



1 2 



A + Zi) + —AaE - aB) 



(22) 



where the mass-matrix V7 j = d'^V / d(pid(pj . The Einstein equations relate the scalar 
metric perturbations to matter perturbations via the energy and momentum con- 
straints (31] 

3H (V> + HA) + + H{a^E - aB)] = -4nG5p , (23) 

iP + HA = -AirGSq , (24) 

where the energy and pressure perturbations and momentum for n scalar fields are 
given by [1] 

Sp = ^ ['Pi (Sfi - 'piA) + Vi5<pi] , (25) 
SP = Y^ [^i {S^i - 'PI A) - ViSifi] , (26) 
Sqa = — ^ (piSifi^i , (27) 

where Vi = dV/dipi. 

We can construct a variety of gauge-invariant combinations of the scalar metric 
perturbations. The longitudinal gauge corresponds to a specific gauge-transformation 
to a (zero-shear) frame such that E — B = 0, leaving the gauge- invariant variables 

$^A~-^^[a\E~B/a)], (28) 
>F = iP + a'^H{E-B/a). (29) 
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Another variable commonly used to describe scalar perturbations during infla- 
tion is the field perturbation in the spatially fiat gauge (where tp — 0). This has the 
gauge-invariant definition [151 116) : 

5(pi^ = 5(fi + ^ip . (30) 

It is possible to use the Einstein equations to eliminate the metric perturbations 
from the perturbed Klein-Gordon equation (|22|l . and write a wave equation solely 
in terms of the field perturbations in the spatially flat gauge |25| 

J 

Only at lowest order in the slow-roll expansion can the interaction terms be 
neglected and we recover the simplifled wave equation (|15|) for a massless field in an 
unperturbed FRW universe. 



Vij 



5^.,^=Q. (31) 



3.3 Adiabatic and entropy perturbations 

There are two more gauge-invariant scalars which are commonly used to describe 
the overall curvature perturbation. The curvature perturbation on uniform-density 
hyp ersurf aces is given by 

-C = V + ^5p, (32) 
P 

first introduced by Bardeen, Steinhardt and Turner [17] (see also Refs. [181 1191 [20] ). 
The comoving curvature perturbation (strictly speaking the curvature perturbation 
on hypersurfaces orthogonal to comoving worldlines) 

n = ^-^^5q, (33) 

where the scalar part of the 3-momentum is given by diSq. TZ has been used by 
Lukash [21[ . Lyth [22[ and many others. For single field inflation we have Sq — —(j)5<j) 
and hence 

n^ii+^Scj). (34) 
The difference between the two curvature perturbations —C, and TZ, 

-(:-TZ= ^Spm , (35) 
P 

is proportional to the comoving density perturbation [23| 

5pm = Sp — 3HSq. (36) 



The energy and momentum constraints ([23p and ([24p can be combined to give a 
generalisation of the Poisson equation 
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relating the longitudinal gauge metric perturbation (|29|l to the comoving density 
perturbation p6[) . Thus the two curvature perturbations, TZ and —(^, coincide on 
large scales {k/aH ^ 1) so long as the longitudinal gauge metric perturbation, <!', 
remains finite - which is generally true during slow-roll infiation. 

The energy conservation equation can be written in terms of the curvature per- 
turbation on uniform-density hypersurfaces, defined in (|32p . to obtain the first-order 
evolution equation |20l [T] 

C = ~H^-S, (38) 
where 5Pnad is the non-adiabatic pressure perturbation, 

<5Pnad = (5P- ^<5p, (39) 
P 

and S is the scalar shear along comoving worldlines [24] , which can be given relative 
to the Hubble rate as 



2p 



9(p + P) a^m 



(40) 



Thus and TZ are constant (and equal upto a sign difference) for adiabatic pertur- 
bations on super-Hubble scales (k/aH <C 1), so long as ^ remains finite, in which 
case the shear of comoving worldlines can be neglected. 

More generally we can define adiabatic perturbations to be perturbations which 
lie along the background trajectory in the phase space of spatially homogeneous fields 
[201 [5] . That is, we generalise Eq. (|39p so that for adiabatic linear perturbations of 
any two variables x and y we require 

Sx Sy 

— = — : adiabatic (41) 
X y 

Thus adiabatic perturbations in a multiple field infiation can be characterised by 
a unique shift along the background trajectory SN = —H5x/x = —H5y/y. For 
example, for adiabatic perturbations of the primordial plasma we require that the 
baryon-photon ratio, UB/n-y, remains unperturbed, and hence 

S{nB/n-,) _ I Shb Sn^, 



3i/p^-:^ =0, (42) 

where we have used hx = —SHrix for baryon number density, and photon number 
density below about 1 MeV. 

For a single scalar field the non-adiabatic pressure (|39|) can be related to the 
comoving density perturbation (|36p [6] 

2V 

5Pnad = -:-77^5pm. (43) 

oil If 



From the Einstein constraint (|37p this will vanish on large scales {k/aH 0) ii >P 
remains finite, and hence single scalar field perturbations become adiabatic in this 
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large-scale limit. In particular, we have from Eq. (|38p that becomes constant for 
adiabatic perturbations in this large scale limit, and hence the curvature perturba- 
tion can be calculated shortly after Hubble-exit in single field inflation and equated 
directly with the primordial curvature perturbation, independently of the details of 
reheating, etc, at the end of inflation. But in the presence of more than one light 
field the vacuum fluctuations stretched to super-Hubble scales will inevitably in- 
clude non-adiabatic perturbations due to the presence of multiple trajectories in the 
phase space. 

We define entropy perturbations to be fluctuations orthogonal to the background 
trajectory 

Sx Su 

Sxy oc : entropy (44) 

X y 

For example in the primordial era we can have entropy perturbations in the primor- 
dial plasma 

S{nB/nj) _ f Sub 5n^ 



Sb = ^ ° = _3Ji- - ^ . (45) 

UB/n^, \nB n-, J 

which are also referred to as baryon isocurvature perturbations. 

In the radiation-dominated era we can define a gauge- invariant primordial curva- 
ture perturbations associated with each of the component fluids [201 126j in analogy 
with the total curvature perturbation (|32|) 

Ci = ~^~ H^4^ , (46) 
Pi 

and these will be constant in the large-scale limit for non-interacting fluids with 
barotropic equation of state Pi{pi) (and hence vanishing non-adiabatic pressure 
perturbations for each fluid) 20,. We can identify isocurvature perturbations, such 
as the baryon isocurvature perturbation (|45|l . with the difference between each C,i 
and, by convention, C,~i for the photons 

Si = 3{(i- G) , (47) 

and the total curvature perturbation ((32]) is given by the weighted sum 



C = E-C- (48) 

Scale-invariant spectra of primordial isocurvature perturbations give rise to dis- 
tinctive power spectra of CMB temperature anisotropies and polarisation |27) and 
are thus tightly constrained by observations |28II29|[30] . although there may be hints 
of isocurvature modes in current CMB data |31) . 

Arbitrary field perturbations in multiple field inflation can be decomposed into 
adiabatic perturbations along the inflaton trajectory and n— 1 entropy perturbations 
orthogonal to the inflaton direction ^ in field space: 

5a = ^ cfi5Lpi , (49) 
&si = ^ sijSipj , (50) 
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where sji&i — 0. Without loss of generality I will assume that the entropy fields 
are also mutually orthogonal in field space. Note that I have assumed that the fields 
have canonical kinetic terms, that is, the field space metric is fiat. See Refs. [7l l32ll33] 
for the generalisation to non-canonical kinetic terms. 

The total momentum and pressure perturbation (|27|) and ((26} for n scalar field 
perturbations can be written in the same form as for a single inflaton field 

Sq — —aSa , (51) 
SP = &{5a - &A) - Va&o . (52) 

However the density perturbation (|25|l is given by 

5p = a(5cr-aA) + K5cr + 25»V, (53) 

where the deviation from the single field result arises due to the non-adiabatic per- 
turbation of the potential orthogonal to the inflaton trajectory: 

bsV = ^ ViSifii - V^So . (54) 
I 

The non-adiabatic pressure perturbation (|39|) is written as [6l[l] 

2V 

SPn.d = -—frSp^. -25,V, (55) 

011(7 

where the comoving density perturbation, Spm, is given by Eq. p6p . Although the 
constraint Eq. (|37p requires the comoving density perturbation to become small on 
large scales, as in the single field case, there is now an additional contribution to the 
non-adiabatic pressure due to non-adiabatic perturbations of the potential which 
need not be small on large scales. 

It is important to emphasise that the presence of entropy perturbations during 
inflation does not mean that the "primordial" density perturbation (at the epoch 
of primordial nucleosynthesis) will contain isocurvature modes. In particular, if the 
universe undergoes conventional reheating at the end of inflation and all particle 
species are driven towards thermal equilibrium with their abundances determined 
by a single temperature (with no non-zero chemical potentials) then the primordial 
perturbations must be adiabatic It is these primordial perturbations that set 
the initial conditions for the evolution of the radiation-matter fluid that determines 
the anisotropics in the cosmic microwave background and large-scale structure in our 
Universe, and thus are directly constrained by observations. We will see that while 
the existence of non-adiabatic perturbations after inflation requires the existence of 
non-adiabatic perturbations during inflation |35j . it is not true that non-adiabatic 
modes during inflation necessarily give primordial isocurvature modes [36j . 



4 Perturbations from two-field inflation 

In this section I will consider the specific example of the coupled evolution of two 
canonical scalar fields, (j) and x, during inflation and how this can give rise to corre- 
lated curvature and entropy perturbations on large scales after inflation [37] . I will 
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use the local rotation in field space defined by Eq. (149 |l and (|50p to describe the 
instantaneous adiabatic and entropy field perturbations. 

The inflaton field perturbation (|49p is gauge-dependent, but we can choose to 
work with the inflaton perturbation in the spatially flat {tfj = 0) gauge: 

Sa^ = 5a + ^tp . (56) 

On the other hand, the orthogonal entropy perturbation (|5Up is automatically gauge- 
invariant. 

The generalisation to two fields of the evolution equation for the inflaton field 
perturbation in the spatially flat gauge, obtained from the perturbed Klein-Gordon 
equations Ip2)l . is [6] 



5<j^ + 3H5a^ + 



, ,^ -2 SttG d / a^&^ 



5a ^ 



and the entropy perturbation obeys 



2±(05s)^2 y^ + ^\ e5s, (57) 



5s + 3H5s +('^ + Vss+ SeA 5s = ^^r^^ , (58) 
\ J a 27rGa^ 

where taxiO = x/0 and 

= (cos' e)V^4, + (sin 20) V^^ + (sin' 6)V-^^, (59) 
Vss = {sin e)V^^ - {sm2e)V^^ + (cos' 9)V^^. (60) 

We can identify a purely adiabatic mode where 5s — on large scales. However a 
non-zero entropy perturbation does appear as a source term in the perturbed inflaton 
equation whenever the inflaton trajectory is curved in field space, i.e., 7^ 0. We 
note that is given by [6] 

0^-^, (61) 
a 

where Vs is the potential gradient orthogonal to the inflaton trajectory in field space. 

The entropy perturbation evolves independently of the curvature perturbation 
on large-scales. It couples to the curvature perturbation only through the gradient 
of the longitudinal gauge metric potential, <!'. Thus entropy perturbations are also 
described as "isocurvature" perturbations on large scales. Eq. (|57p shows that the 
entropy perturbation 5s works as a source term for the adiabatic perturbation. This 
is in fact clearly seen if we take the time derivative of the curvature perturbation 

m-. 

TZ ^iL'^^+^eSs, (62) 
Ha^ a 

Therefore TZ (or ^) is not conserved even in the large-scale limit in the presence of 
the entropy perturbation 5s with a non-straight trajectory in field space (0^0). 

Analogous to the single field case we can introduce slow-roll parameters for light, 
weakly coupled fields [38]. At first-order in a slow-roll expansion, the inflaton rolls 
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directly down the potential slope, that is Vs — 0. Thus we have only one slope 
parameter 

but three parameters, rjaa, rjas and rjss, describing the curvature of the potential, 
where 

The background slow-roll solution is described in terms of the slow-roll param- 
eters by 

a^-^el/, H'^c^^n^s, (65) 
while the perturbations obey 

H'^Sa^ ~ (2e — rja-a-) Sa^ — 2Tjaa5s , 

H'^Ss ~ —rjssSs , (66) 

on large scales, where we neglect spatial gradients. Although 14 — at lowest order 
in slow-roll, this does not mean that the inflaton and entropy perturbations decouple. 
6 given by Eq. (|65p is in general non-zero at first-order in slow-roll and large-scale 
entropy perturbations do affect the evolution of the adiabatic perturbations when 

??<TS ^ 0. 

While the general solution to the two second-order perturbation equations (|57[l 
and (|58l) has four independent modes, the two first-order slow-roll equations (|66[l 
give the approximate form of the squeezed state on large scales. This has only two 
modes which we can describe in terms of dimensionless curvature and isocurvature 
perturbations: 

U H 
TZ = — 5a^ , S= — 5s. (67) 
a a 

The normalisation of TZ coincides with the standard definition of the comoving curva- 
ture perturbation, Eq. (|33p . The normalisation of the dimensionless entropy during 
inflation, S, is chosen here coincide with Ref. fSS]. It can be related to the non- 
adiabatic pressure perturbation (I39|l on large scales 

5Pnad ^ -^V<rs-;—^S . (68) 

zttG 

The slow-roll approximation can provide a useful approximation to the instanta- 
neous evolution of the flelds and their perturbations on large scales during slow-roll 
inflation, but is not expected to remain accurate when integrated over many Hubble 
times, where inaccuracies can accumulate. In single-field inflation the constancy of 
the comoving curvature perturbation after Hubble exit, which does not rely on the 
slow-roll approximation, is crucial in order to make accurate predictions of the pri- 
mordial perturbations using the slow-roll approximation only around Hubble cross- 
ing. In a two-fleld model we must describe the evolution after Hubble exit in terms 
of a general transfer matrix: 



TZ\ ^ f I Tns\ f TZ 
S j \QTss )\S 



(69) 



On large scales the comoving curvature perturbation still remains constant for the 
purely adiabatic mode, corresponding to 5 = 0, and adiabatic perturbations remain 
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adiabatic. These general results are enough to fix two of the coefficients in the 
transfer matrix, but T-jzs and Tss remain to be determined either within a given 
theoretical model, or from observations, or ideally by both. The scale-dependence 
of the transfer functions depends upon the inflaton-entropy coupling at Hubble exit 
during inflation and can be given in terms of the slow-roll parameters as [38] 

d 

T-R.S = 2r;CTs + (2e — 7]^^ + riss)Tns , 

Tss = (2e - + Vss)Tss ■ (70) 



dink 

d 

d\nk 



4.1 Initial power spectra 

For weakly-coupled, light flelds we can neglect interactions on wavelengths below the 
Hubble scale, so that vacuum fluctuations give rise to a spectrum of uncorrelated 
field fluctuations on the Hubble scale {k = aH) during inflation given by Eq. (|17p : 

(^)' , (71) 

where we use a * to denote quantities evaluated at Hubble-exit. If a field has a mass 
comparable to the Hubble scale or larger then the vacuum fiuctuations on wave- 
lengths greater than the effective Compton wavelength are suppressed. In addition 
fiuctuations in strongly interacting fields may develop correlations before Hubble 
exit. But during slow-roll infiation the correlation between vacuum fiuctuations in 
weakly coupled, light fields at Hubble-exit is suppressed by slow-roll parameters. 
This remains true under a local rotation in fields space to another orthogonal ba- 
sis such as the instantaneous infiaton and entropy directions (|49|) and (|50p in field 
space. 

The curvature and isocurvature power spectra at Hubble-exit are given by 
while the cross-correlation is first-order in slow-roll |49l 150) . 

CTZs\,~-2Cna.VTl\,, (73) 

where C = 2 — ln2 — 7«0.73 and 7 is the Euler number. The normalisation chosen 
for the dimensionless entropy perturbation in Eq. (|67p ensures that the curvature 
and isocurvature fiuctuations have the same power at horizon exit [38]. The spectral 
tilts at horizon-exit are also the same and are given by 

AnnU^AnsU~-6e + 2r]^^. (74) 

where Anx = dlnVx / d\nk. 

The tensor perturbations (|2U|I are decoupled from scalar metric perturbations 
at first-order and hence the power spectrum has the same form as in single field 
infiation. Thus the power spectrum of gravitational waves on super-Hubble scales 
during infiation is given by 
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_ im' ^ 128 v; 

and the spectral tilt is 

Am:\,~-2e. (76) 



4.2 Primordial power spectra 

The resulting primordial power spectra on large scales can be obtained simply by 
applying the general transfer matrix (|69p to the initial scalar perturbations. The 
scalar power spectra probed by astronomical observations are thus given by 08] 

VTz = {l+Tls)VnU (77) 
Vs = TlsVnU (78) 
Ctis = TtzsTssVtz]* . (79) 

The cross-correlation can be given in terms of a dimensionless correlation angle: 

rn.P,= ^Jl^= (80) 

We see that if we can determine the dimensionless correlation angle, 0, from 
observations, then this determines the off-diagonal term in the transfer matrix 

Ttis = cot © , (81) 

and we can in effect measure the contribution of the entropy perturbation during 
two-field infiation to the resultant curvature primordial perturbation. In particular 
this allows us in principle to deduce from observations the power spectrum of the 
curvature perturbation at Hubble-exit during two-field slow-roll inflation |38) : 

VTz\*^Vns\n'e. (82) 

The scale-dependence of the resulting scalar power spectra depends both upon 
the scale-dependence of the initial power spectra and of the transfer coefficients. 
The spectral tilts are given from Eqs. (|77H79|) by 

Ann = AnnU + H'^ {dTns / dU) sin 20 , 

Ans = AnnU + 2H:\d\nTss / dU) , (83) 
Anc = AmzU + [{dTns/dU)ta.n0 + (dlnTss/dU)] , 

where we have used Eq. (|8ip to eliminate T-jis in favour of the observable correlation 
angle 0. Substituting Eq. (|74|) for the tilt at Hubble-exit, and Eqs. (|70[) for the scale- 
dependence of the transfer functions, we obtain [3S] 

Arm ^ -(6- 4cos^6')e 

+2 (r]a-cr sin^ + 2r]a-s sin cos + rjgs cos^ , 
Ans ^ -2e + 2tJss , (84) 
Anc ~ — 2e -I- 2riss + 2r]as tan . 
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Although the overall amplitude of the transfer functions are dependent upon 
the evolution after Hubble-exit and through reheating into the radiation era, the 
spectral tilts can be expressed solely in terms of the slow-roll parameters at Hubble- 
exit during inflation and the correlation angle, 0, which can in principle be observed. 

If the primordial curvature perturbation results solely from the adiabatic inflaton 
field fluctuations during inflation then we have T-jzs = in Eq. (|77p and hence 
cos© = in Eqs. (|84p . which yields the standard single field result 

Ann ^ -6e + 2rj^„ . (85) 

Any residual isocurvature perturbations must be uncorrelated with the adiabatic 
curvature perturbation (at first-order in slow-roll) with spectral index 

Ans ~ -2e + 2riss ■ (86) 

On the other hand, if the observed primordial curvature perturbation is produced 
due to some entropy field fluctuations during inflation, we have T-ns 3> 1 and sin ~ 
0. In a two-field inflation model any residual primordial isocurvature perturbations 
will then be completely correlated (or anti-correlated) with the primordial curvature 
perturbation and we have 

Ann ^ Anc ~ Ans ~ -2e + 2i]ss ■ (87) 

The gravitational wave power spectrum is frozen-in on large scales, independent 
of the scalar perturbations, and hence 

Vt=Vt\,. (88) 

Thus we can derive a modified consistency relation [39] between observables appli- 
cable in the case of two-field slow-roll infiation: 

r = ^ ~ -8AnT sin . (89) 

This relation was first obtained in Ref. [40) at the end of two-field infiation, and 
verified in Ref. [41) for slow-roll models. But it was realised in Ref. [38] that this 
relation also applies to the primordial perturbation spectra in the radiation era long 
after two-field slow-roll infiation has ended and hence may be tested observationally. 

More generally, if there is any additional source of the scalar curvature per- 
turbation, such as additional scalar fields during inflation, then this could give an 
additional contribution to the primordial scalar curvature spectrum without affect- 
ing the gravitational waves, and hence the more general result is the inequality [42| : 

r < -8AnT sin^ . (90) 

This leads to a fundamental difference when interpretting the observational con- 
straints on the amplitude of primordial tensor perturbations in multiple inflation 
models. In single field infiation, observations directly constrain r = [Pt/Vtz]* and 
hence, from Eqs. ()72)) and ()75[l . the slow-roll parameter e. However in multiple field 
inflation, non-adiabatic perturbations can enhance the power of scalar perturbations 
after Hubble exit and hence observational constraints on the amplitude of primordial 
tensor perturbations do not directly constrain the slow-roll parameter e. 

Current CMB data alone require r < 0.55 (assuming power-law primordial spec- 
tra) [43[ which in single-field models is interpretted as requiring e < 0.04. But in 
multiple field models e could be larger if the primordial density perturbation comes 
from non-adiabatic perturbations during inflation. 
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5 Non-Gaussianity 

A powerful technique to calculate the primordial curvature perturbation resulting 
from many inflation models, including multi-fleld models, is to note that the curva- 
ture perturbation ( defined in Eq. (|32p can be interpreted as a perturbation in the 
local expansion [44] |42] |46] 

C = SN, (91) 

where 5N is the perturbed expansion to uniform-density hypersurfaces with respect 
to spatially flat hypersurfaces, which is given to first-order by 

^ = , (92) 

P 

where 5p^ must be evaluated on spatially flat {tp — 0) hypersurfaces. 

An important simpliflcation arises on large scales where anisotropy and spatial 
gradients can be neglected, and the local density, expansion, etc, obeys the same 
evolution equations as in a homogeneous FRW universe [421 1451 1201 1241 1471 146] . 
Thus we can use the homogeneous FRW solutions to describe the local evolution, 
which is known as the "separate universe" approach [421 1451 1201 147| . In particular 
we can evaluate the perturbed expansion in different parts of the universe resulting 
from different initial values for the fields during inflation using the homogeneous 
background solutions [42]. The integrated expansion from some initial spatially fiat 
hypersurface up to a late-time fixed density hypersurface, say at the epoch of primor- 
dial nucleosynthesis, is some function of the field values on the initial hypersurface, 
N{ipi\^). The resulting primordial curvature perturbation on the uniform-density 
hypersurface is then 

I 

where Nj = dN/dipi and Stpj^ is the field perturbation on some initial spatially- 
flat hypersurfaces during inflation (|30p . In particular the power spectrum for the 
primordial density perturbation in a multi-fleld inflation can be written (at leading 
order) in terms of the held perturbations after Hubble-exit as 

n^^iVlP,^,^. (94) 

This approach is readily extended to estimate the non-linear effect of field per- 
turbations on the metric perturbations [451 1241 146] . We can take Eq. ([91[) as our 
definition of the non-linear primordial curvature perturbation, (, so that in the ra- 
diation dominated era the non-linear extension of Eq. ([92p is given by [46] 



C = 7ln(-) , (95) 



4 yp 

where p{t, x) is the perturbed (inhomogeneous) density evaluated on a spatially fiat 
hypersurface and p{t) is the background (homogeneous) density. This non-linear 
curvature perturbation as a function of the initial field fiuctuations can simply be 
expanded as a Taylor expansion [571 1551 1481 152[ 



C — ^ NiSifi^ + i ^ NijSipi^S<pj^ + i ^ NijKS(pi^S<fij^5<fiKi, + . . . . (96) 



2^ "■'^'^^•"^ , g 

I.J I,J,K 



18 David Wands 



where we now identify (|93p as the leading-order term. 

We expect the field perturbations at Hubble-exit to be close to Gaussian for 
weakly coupled scalar fields during infiation [531 1541 1561 1551 151| . In this case the 
bispectrum of the primordial curvature perturbation at leading (fourth) order, can 
be written using the iJA'^-formalism, as [751 157] 



Bc(ki,k2,k3) = pNL[Pdki)Pc{k2) + Pdk2)Pc{k:i) + Pdk3)Pc{ki)] ■ (97) 

where -P^(fc) = 2n^P^{k)/k^ , and the dimensionless non-linearity parameter is given 
by [57] 

^ 6 (NcNCf ^ ' 

Similarly, the connected part of the trispectrum in this case can be written as [521148] 

rc(ki,k2,k3,k4) = TNL [Pc(|ki +k3|)Pc(fc3)Pc(fc4) + (U pemis)] 
54 

+ ^9NL [Pc{k2)Pc{k3)Pdk4) + (3 perms)] . (99) 

where 

= (jv,,jv^)3 ■ (^00) 

25NabcN'^N''N'=' .^^^^ 
- 54 [NnN^r ■ ^ ^ 

The expression for tml was first given in [58] . Note that we have factored out 
products in the trispectrum with different k dependence in order to define the two k 
independent non-linearity parameters rjvz, and Qnl- This gives the possibility that 
observations may be able to distinguish between the two parameters [59] . 

In many cases there is single direction in field-space, x, which is responsible for 
perturbing the local expansion, N{x), and hence generating the primordial curvature 
perturbation ([96]1 . For example this would be the inflaton field in single field models 
of inflation, or it could be the late-decaying scalar field in the curvaton scenario 
]61l 1621 163] as will be discussed in the next section. In this case the curvature 
perturbation (]96p is given by 

C ~ N'Sxi, + ^N"Sxl + In"'5x^p + ..., (102) 

and the non-Gaussianity of the primordial perturbation has the simplest "local" 
form 

C = Ci + |/iVLC? + ;^3iVLC? + ... (103) 
5 25 

where C^i = N'Sx^i is the leading-order Gaussian curvature perturbation and the 
non-linearity parameters /jvl and gNL, are given by [571 [M] 



5 iV" 



25 A^'" 

^^^ = 54(iV03' (1°^) 
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The primordial bispectrum and trispectrum are then given by Eqs. (|97|l and (|99|l . 
where the non-linearity parameters /jvz, and gNL, given in Eqs. (|98|l and (|101|l . 
reduce to Eqs. (I104|l and (|105|) respectively, and tnl given in Ea. (|100|) reduces to 

Thus TjvL is proportional to /^^ (first shown in [59) using the Bardeen potential, 
and in [57] using this notation) . However the trispectrum could be large even when 
the bispectrum is small because of the qnl term [591 164] . 

In the case of where the primordial curvature perturbation is generated solely 
by adiabatic fluctuations in the inflaton fleld, a, the curvature perturbation is non- 
linearly conserved on large scales [761 1461 [50] and we can calculate A'^', A'^", A'^'", etc, 
at Hubble-exit. In terms of the slow-roll parameters, we find 

N' = ^L^'±'^0(e-^), (107) 

N" ^~^-^-(v..-2e)^0{l) , (108) 

N'" ^^^-^ (eiiaa - ilL + Ul) ~ 0{e^) , (109) 

where we have used the reduced Planck mass = (SttG)"^ and introduced the 
second-order slow-roll parameter — ra^VaVacrcr jV'^ ■ Hence the non-linearity pa- 
rameters for single field infiation, (I104[l and ([105p . are given by 



/iVL = |(r?<..-2e), (110) 
gNL = ^ (2e77CT<T - 2ril„ + . (Ill) 

with Tnl given by Eq. ([106p . Although there are additional contributions to the 
primordial bispectrum and trispectrum coming from the intrinsic non-Gaussianity 
of the field perturbations at Hubble-exit, these are also suppressed by slow-roll 
parameters in slow-roll inflation. Thus the primordial non-Gaussianity is likely to 
be too small to ever be observed in the conventional inflaton scenario of single-fleld 
slow-roll inflation. Indeed any detection of primordial non-Gaussianity fNL > 1 
would appear to rule out this inflaton scenario. 

However significant non-Gaussianity can be generated due to non-adiabatic field 
fiuctuations. Thus far it has proved difficult to generate significant non-Gaussianity 
in the curvature perturbation during slow-roll infiation, even in multiple field models. 
But detectable non-Gaussianity can be produced when the curvature perturbation 
is generated from isocurvature field perturbations at the end of infiation [651 112] , 
during inhomogeneous reheating [66II67IIT4] . or after infiation in the curvaton model, 
which I will discuss next. 



6 Curvaton scenario 



Consider a light, weakly-coupled scalar field, X) that decays some time after infiation 
has ended. There are many such scalar degrees of freedom in supersymmetric theories 
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and if they are too weakly coupled, and their lifetime is too long, this leads to the 
"Polonyi problem" [71) . Assuming the field is displaced from the minimum of its 
effective potential at the end of inflation, the field evolves little until the Hubble rate 
drops below its effective mass. Then it oscillates, with a time-averaged equation of 
state of a pressureless fluid, = 0, (or, equivalently, a collection of non-relativistic 
particles) and will eventually come to dominate the energy density of the Universe. 
To avoid disrupting the standard "hot big bang" model and in particular to preserve 
the successful radiation-dominated model of primordial nucleosythesis, we require 
that such fields decay into radiation before t ~ 1 second. For a weakly-coupled 
field that decays with only gravitational strength, F ~ m'^/Mp, this requires > 
100 TeV. 

But there is a further important feature of late-decaying scalar fields that has 
only recently received serious consideration. If the field is inhomogeneous then it 
could lead to an inhomogeneous radiation density after it decays [721 173] . This is 
the basis of the curvaton scenario [611 1621 [63] . 

If this field is light (m < H) during inflation then small-scale quantum fluctu- 
ations will lead to a spectrum of large-scale perturbations, whose initial amplitude 
at Hubble-exit is given (at leading order in slow-roll) by Eq. ([17[) . When the Hubble 
rate drops and the field begins oscillating after inflation, this leads to a flrst-order 
density perturbation in the x-field: 

Cx = -V'+|^. (112) 

where = Cx remains constant for the oscillating curvaton field on large 

scales, so long as we can neglect its energy loss due to decay. Using Eq. I|17p for the 
field fluctuations at Hubble-exit and neglecting any non-linear evolution of the x- 
field after inflation (consistent with our assumption that the field is weakly coupled) , 
we have 

The total density perturbation (|32|l . considering radiation, 7, and the curvaton, x, 
is given by [62] 

^ 4p-,Cj + 3pyCx 
4p^+3px 

Thus if the radiation generated by the decay of the inflaton at the end of inflation is 
unperturbed {V^''^ <C 10~^) the total curvature perturbation grows as the density 
of the x-fi^ld grows relative to the radiation: ~ J^xCx- 

Ultimately the x-field must decay (when H ^ F) and transfer its energy den- 
sity and, crucially, its perturbation to the radiation and/or other matter fields. In 
the simplest case that the non-relativistic x-field decays directly to radiation a full 
analysis [701 174[ of the coupled evolution equation gives the primordial radiation 
perturbation (after the decay) 

C = r^{p)(,x , (115) 
where p = [^2x/(-^/^)^^^]initiai is a dimensionless parameter which determines the 

m 

n Q24 \ -1-24 

,(p) ^ 1 - ( 1 + . (116) 



Multiple field inflation 21 



For p 3> 1 the x-field dominates the total energy density before it decays and r^^ ~ 1, 
while for p <C 1 we have ~ 0.924p <^ 1. 

Finally combining Eqs. (|113|l and (|115|l we have 



Note that the primordial curvature perturbation tends to have less power on small 
scales due to the decreasing Hubble rate at Hubble exit in Eq. (|117|) . but can also 
have more power on small scales due to the decreasing for a positive efi'ective 
mass-squared, during inflation. In terms of slow-roll parameters the actual tilt is 
given by Eq. (|84|l when sin — 



In the extreme slow-roll limit the spectrum becomes scale- invariant, as in the inflaton 
scenario. 

In contrast to the inflaton scenario the flnal density perturbation in the curvaton 
scenario is a very much dependent upon the physics after the field perturbation 
exited the Hubble scale during infiation. For instance, if the curvaton lifetime is too 
short then it will decay before it can significantly perturb the total energy density 
and V^^^ <^ 10~^. The observational constraint on the amplitude of the primordial 
perturbations gives a single constraint upon both the initial fluctuations during 
inflation and the post-inflationary decay time. This is in contrast to the inflaton 
scenario where the primordial perturbation gives a direct window onto the dynamics 
of inflation, independently of the physics at lower energies. In the curvaton scenario 
there is the possibility of connecting the generation of primordial perturbations to 
other aspects of cosmological physics. For instance, it may be possible to identify 
the curvaton with flelds whose late-decay is responsible for the origin of the baryon 
asymmetry in the universe, in particular with sneutrino models of leptogenesis (in 
which an initial lepton asymmetry is converted into a baryon asymmetry at the 
electroweak transition) [68]. 

The curvaton scenario has re-invigorated attempts to embed models of inflation 
in the very early universe within minimal supersymmetric models of particle physics 
constrained by experiment It may be possible that the inflaton field driving 
inflation can be completely decoupled from visible matter if the dominant radiation 
in the universe today comes from the curvaton decay rather than reheating at the 
end of inflation. Indeed the universe need not be radiation-dominated at all until 
the curvaton decays if instead the inflaton fast-rolls at the end of inflation. 

The curvaton offers a new range of theoretical possibilities, but ultimately we 
will require observational and/or experimental predictions to decide whether the 
curvaton, inflaton or some othe field generated the primordial perturbation. 1 will 
discuss observational predictions of the curvaton scenario in the following subsection. 

6.1 Non-Gaussianity 

The best way to distinguish between different scenarios for the origin of structure 
could be the statistical properties of the primordial density perturbation. Primordial 
density perturbations in the curvaton scenario originate from the small-scale vacuum 




(117) 



Ann ^ -2e + 2ri- 



'XX ■ 



(118) 
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fluctuations of the weakly interacting curvaton field during inflation, which can 
be described on super-Hubble scales by a Gaussian random field. Thus deviations 
from Gaussianity in the primordial bispectrum and connected trispectrum can be 
parameterised by the dimensionless parameters /jvl and qnl defined in Eqs. (|104p 
and (fT05l) . 

When the curvaton field begins oscillating about a quadratic minimum of its 
potential we have Px = and the time-averaged equation of state becomes, 



Px = 0. The non-linear generalisation of the primordial curvature perturbation (|46|) 
on hypersurfaces of uniform-curvaton density is then 

where we distinguish here between the inhomogeneous density on spatially fiat 
hypersurfaces and the average density p-^- Given that cx we thus have 

Cx^mfi + ^^^l . (120) 



This gives the full probability distribution function for C^^ for Gaussian field pertur- 
bations 5x- Expanding to first-order we obtain C|"xi = 2Sx-4:/3x, and then to second- 
and third-order we obtain by analogy with Eq. (I103|l the non-linearity parameters 
for the curvaton density perturbation 

Cx ^ Cxi + f /^lCxI + ^9^l(x1 + • ■ • . (121) 

where [64] 

f^L - -I , (122) 

91l-§- (123) 

If the curvaton dominates the total energy density before it decays into radiation, 
then this is the curvature perturbation, and specifically the non-linearity parameters, 
inherited by the primordial radiation density. Although not suppressed by slow-roll 
parameters, this non-Gaussianity is still smaller than the best upper limits expected 
from the Planck satellite [75j . 

On the other hand if the curvaton decays before it dominates over the energy 
density of the existing radiation, so the transfer function r^ip) <C 1 in Eq. (|115p . then 
the curvaton may lead to a large and detectable non-Gaussianity in the radiation 
density after it decays. Assuming the sudden decay of the curvaton on the H — 
r uniform-density hypersurface leads to a non-linear relation between the local 
curvaton density and the radiation density before and after the decay [64) 

p^e-^^ + Pxe^^'^--^^ =p^ + Px. (124) 
Expanding this term by term yields [771 1571 164) 

^--4^-f-T^' (^^^) 

... = -^ + ^ + i^ + ^. (126) 
^ 6ry 108 27 18 ^ ^ 
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These reduce to Eqs. (|122|l and (|123|l as r;^ ^ 1, but become large for <C 1. 
These analytic results rely on the sudden decay approximation but have been tested 
against numerical solutions [781 164) and give an excellent approximation for both 
<C 1 and r-^ ~ 1. 

More generally one can use Eqs. (|119p and (|124p . or use Eq. (|91|l and solve the 
non-linear, but homogeneous equations of motion to determine N(x) to give the full 
probability distribution function for the primordial curvature perturbation ^ |64j . 

Current bounds from the WMAP sateUite require — 54 < Jml < 114 at the 95% 
confidence limit [43], and hence require > 0.011. But future cosmic microwave 
background (CMB) experiments such as Planck could detect Jnl as small as around 
5 [75] , and it has even been suggested that it might one day be possible to constrain 
/jvi ~ 0.01 [79]. 

6.2 Residual isocurvature perturbations 

In the curvaton scenario the initial curvaton perturbation is a non-adiabatic pertur- 
bation and hence can in principle leave behind a residual non-adiabatic component. 
Perturbations in this one field, would be responsible for both the total primordial 
density perturbation and any isocurvature mode and hence there is the clear pre- 
diction that the two should be completely correlated, corresponding to cos© = ±1 
in Eq. JHO} and Arm - 1 = Ans = Anc in Eqs. (|84|l . 

Using C,i defined in Eq. (|46p for different matter components it is easy to see 
how the curvaton could leave residual isocurvature perturbations after the curvaton 
decays. If any fluid has decoupled before the curvaton contributes significantly to 
the total energy density that fiuid remains unperturbed with C,i — 0, whereas after 
the curvaton decays into radiation the photons perturbation is given by (|115|) . Thus 
a residual isocurvature perturbation (|47p is left 

5i = -3C, (127) 

which remains constant for decoupled perfect fluids on large scales. 

The observational bound on isocurvature matter perturbations completely cor- 
related with the photon perturbation, is [80] 

- 0.42 < + (Pcdm/pB)g,dm ^ ^^28) 

In particular if the baryon asymmetry is generated while the total density perturba- 
tion is still negligible then the residual baryon isocurvature perturbation, Sb = — 3(^.y 
would be much larger than the observational bound and such models are thus ruled 
out. The observational bound on CDM isocurvature perturbations are stronger by 
the factor pcdm/ps .81^ although CDM is usually assumed to decouple relatively 
late. 

An interesting amplitude of residual isocurvature perturbations might be realised 
if the decay of the curvaton itself is the non-equilibrium event that generates the 
baryon asymmetry. In this case the net baryon number density directly inherits the 
perturbation C^b = Cx while the photon perturbation C^-y < may be diluted by 
pre-existing radiation and is given by Eq. (|115p . Note that so long as the net baryon 
number is locally conserved it defines a conserved perturbation on large scales. 
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even though it may still be interacting with other fluids and fields [75]. Hence the 
primordial baryon isocurvature perturbation (|45|) in this case is given by 

5fl=3(l-rJCx = ^^i^^C7- (129) 

Thus the observational bound (|128p requires r-^ > 0.92 if the baryon asymmetry is 
generated by the curvaton decay. 

There is no lower bound on the predicted amplitude of residual non-adiabatic 
modes and, although the detection of completely correlated isocurvature perturba- 
tions would give strong support to the curvaton scenario, the non-detection of pri- 
mordial isocurvature density perturbations cannot be used to rule out the curvaton 
scenario. In particular, if the curvaton decays at sufficiently high temperature and all 
the particles produced relax to a thermal equilibrium abundance, characterised by 
a common temperature (and vanishing chemical potentials) then no residual isocur- 
vature perturbations survive. In full thermal equilibrium there is a unique attractor 
trajectory in phase-space and only adiabatic perturbations (along this trajectory) 
survive on large scales. 



7 Conclusions 

Infiation offers a beautifully simple origin for structure in our Universe. The zero- 
point fluctuations of the quantum vacuum state on sub-atomic scales are swept up 
by the accelerated expansion to astronomical scales, seeding an almost Gaussian 
distribution of primordial density perturbations. The large scale structure of our 
Universe can then form simply due to the gravitational instability of overdense 
regions. 

Astronomical observations over recent years have given strong support to this 
simple picture. But increasingly precise astronomical data will increasingly allow 
us to probe not only the parameters of what has become the standard cosmologi- 
cal model, but also to probe the nature of the primordial perturbations from which 
structure formed. Any evidence of primordial gravitational waves, primordial isocur- 
vature fluctuations, and/or non-Gaussianity of the primordial perturbations could 
provide valuable information about the inflationary dynamics that preceded the hot 
big bang. 

Single-fleld slow-roll inflation predicts adiabatic density perturbations with neg- 
ligible non-Gaussianity, but could produce a gravitational wave background which 
could be detected by upcoming CMB experiments. 

On the other hand multiple fleld inflation can lead to a wider range of possi- 
bilities which could be distinguished by observations. A spectrum of non-adiabatic 
field fluctuations on large scales during inflation could leave residual isocurvature 
perturbations after inflation, which can be correlated with the primordial curva- 
ture perturbation, and can give rise to a detectable level of non-Gaussianity. In 
simple models, such as the curvaton scenario, where the primordial curvature per- 
turbations originate from almost Gaussian fluctuations in a single scalar fleld, any 
residual isocurvature perturbations are expected to be completely correlated with 
the curvature perturbation and the non-Gaussianity is of a speciflc "local" form. 
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After 25 years studying inflation, we may for the first time have evidence of 
a weak scale-dependence of the power spectrum of the primordial curvature per- 
turbation [43] which would begin to reveal the slow-roll dynamics during inflation. 
Primordial perturbations may have much more to tell us about the physics of infla- 
tion in the future. 
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